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Abstract. We use large scale Monte Carlo simulations to test scaling theories for the electrostatic persis- 
tence length l e of isolated, uniformly charged polymers with Debye-Hiickel intrachain interactions in the 
limit where the screening length k _1 exceeds the intrinsic persistence length of the chains. Our simulations 
cover a significantly larger part of the parameter space than previous studies. We observe no significant 



deviations from the prediction l e 



by Khokhlov and Khachaturian which is based on applying the 



Odijk-Skolnick-Fixman theory to the stretched de Gennes-Pincus-Velasco-Brochard polyelectrolyte blob 
chain. A linear or sublinear dependence of the persistence length on the screening length can be ruled out. 
We argue that previous numerical results pointing into this direction are probably due to a combination 
of excluded volume and finite chain length effects. The paper emphasizes the role of scaling arguments in 
the development of useful representations for experimental and simulation data. 
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1 Introduction 



The theoretical understanding of macromolecules carry- 
' ing ionizable groups is far from complete 0,0. In spite of 
, the long range of the interactions, the systems are often 
discussed using analogies to neutral polymers. A promi- 
nent example is the concept of an electrostatic persis- 
tence length, which was introduced more than 20 years 
■ ago by Odijk [| and by Skolnick and Fixman § (OSF). 
They considered a semiflexible polymer or wormlike chain 
(WLC) with intrinsic persistence length Iq and Debye- 
Hiickel (DH) screened electrostatic interactions Udh/^bT - 
(q 2 lB/r) exp(-nr) between charges q/e spaced at regular 
intervals A along the chain. The Bjcrrum length lg char- 
acterizes the strength of the electrostatic interactions and 
• is defined as the distance where the Coulomb energy of 
two unit charges e is equal to kgT. Due to the presence 
of mobile ions the bare Coulomb interaction is cut off be- 
yond the screening length kT 1 . OSF were interested in 
bending fluctuations and considered the resulting increase 
of the electrostatic energy relative to the straight ground 
state. In the ("OSF") limit k" 1 <C Iq where the screening 
length is smaller than the intrinsic persistence length of 
the chain and to lowest order in the local curvature, the 
Debye-Hiickel interaction makes an additive contribution 
to the bending rigidity. As a consequence, a WLC with 



DH interactions (DHWLC) behaves in this limit on large 
length scales like an ordinary WLC with renormalized per- 
sistence length 



lp — Iq 
loSF = 



f loSF 
q 2 l B 

AA 2 k 2 



(1) 
(2) 
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Ever since, there has been a lively debate on how to ex- 
tend the theory to parameter ranges beyond the OSF limit 
k _1 -C Iq- Barrat and Joanny (BJ)|| have shown that the 
' original OSF derivation breaks down, if the chains start 
to bend significantly on length scales comparable to the 
screening length. As a consequence, Eq. (|2|) cannot sim- 
ply remain valid beyond the OSF limit as was sometimes 
speculated Two main scenarios, which we denote by 
"OSFKK " and "KMBJ " after the initials of the main 
authors, have been discussed in the literature: 

OSFKK According to Khokhlov and Khachaturian (KK) jj 
the OSF theory can be applied to a "stretched chain 
of polyelectrolyte blobs" , a concept introduced by de 
Gennes et al. (8) to describe the behavior of weakly 
charged flexible polyelectrolytes in the absence of screen- 
ing. The persistence length of the blob chain is then 
calculated from Eq. (j^) using suitably renormalized 
parameters. 

KMBJ Refinements by Muthukumar @[H|[ll| °f the orig- 
inal theory of Katchalsky |lj| treat electrostatic inter- 
actions in strict analogy to short-range excluded vol- 
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ume interactions. Quite interestingly, the results are 
consistent with the scaling picture of de Gennes et 
al. H in the two limits of strong and vanishing screen- 
ing. Moreover, they are supported by recent calcula- 
tions by BJ and others [6, 
sistence length of the bio' 
cedure and found L ~ kT 



13 who determined the per- 



d chain in a variational pro- 



In addition, there is a number of recent theories which fall 
into neither of the two classes outlined above [ ^4| , p^ ,^6|,p7|, 
|Ls| ]. While there is a growing consensus among theoreti- 
cians that the OSFKK result is asymptotically correct jl9[ 
pfj|,^l[, experiment s [p2[ p^,pj|,p5|,p6| as well as computer 
simulations 121 



,|28|,|29j have consistently provided evidence 
for a comparatively weak K-dependence of the electrostatic 
presistence length. 

The purpose of the present paper is to shed some new 
light on this problem by combining a scaling analysis with 
large scale Monte Carlo simulations. We reexamine the 
KMBJ and the OSFKK theory in order to extract guid- 
ance for the data analysis and the choice of simulation 
parameters. As a result we are able (i) to rule out the 
KMBJ theory, (ii) to provide benchmark results for ana- 
lytical solutions of the DHWLC model as well as (iii) for 
a comparison to experiments in order to clarify if the sys- 
tems under consideration are actually described by the 
DHWLC model or if additional effects such as solvent 
quality or counter-ion condensation need to be taken into 
account as well. 

The paper is organized as follows: In section |^ we re- 
view the predictions of the KMBJ and OSFKK theories, 
followed by a discussion in section || of how experiments 
and simulations should be set up and analyzed in order to 
discriminate between the two scaling pictures. Details of 
our Monte Carlo simulations can be found in in section [|. 
We present our results in section || and close with a brief 
discussion. 



2 Scaling theories of intrinsically flexible 
Debye-Hiickel chains 

The DHWLC is characterized by the following set of pa- 
rameters: q, Ib, k, A, /, and L to t, where L to t denotes 
the total chain length and / < 1 the fraction of ionized 
charged groups which needs to be determined indepen- 
dently for experimental systems. Setting / = Ib/A is a 
rough way of accounting for Manning condensation jl]] of 
counter ions in cases where Ib < A. Correlation func- 
tions can be calculated for chain segments with arbitrary 
contour length L < L to t- We focus on the non-OSF limit 
Ib,A, l <C <C \J (r 2 ) (L t ot) where the screening length 
is larger than all microscopic length scales of the polymer 
model but smaller than the size of the entire chain. 

For an understanding of the physics, some of these 
length scales and parameters are less relevant than others. 
For example, the actual distribution of the charges on the 
chain should be unimportant as long as A <C k . In our 
simulations we use discrete charges spaced by a distance 
equal to the intrinsic persistence length, while the scaling 



arguments assume a continuous charge distribution. Sim- 
ilarly, in the non-OSF limit with l <C n^ 1 the physics 
should not depend on the details of the WLC crossover 
from rigid rod to random coil behavior for L w Iq- In our 
simulations we therefore use freely jointed chains (FJC) 
whose (Kuhn) bond length b corresponds (up to a hence- 
forth neglected factor of two) to the persistence length of a 
WLC. Finally, it is convenient to consider the limit of infi- 
nite total chain length in order to eliminate the L to t depen- 
dence. Again the practical limitation to N to t = L to t/b = 
4096 segments should be unimportant, since our chains 
always fulfill the condition k^ 1 J (r 2 ) (L to t)- 

The remaining independent parameters (the line charge 
density fq/A, Ib, k, Iq and L) can be reduced further us- 
ing the notion of a "polyelectrolyte blob" which was in- 
troduced by de Gennes et al. ||] to describe the crossover 
from locally unchanged chain statistics to stretching on 
long length scales. 

Consider first weakly charged flexible polyelectrolytes, 
where the electrostatic interactions are irrelevant on the 
length scales comparable to the intrinsic persistence length 
lo. On larger length scales an undisturbed WLC with a 
contour length L has a spatial extension (r 2 ^ — 2LIq. 
Neglecting prefactors, the electrostatic energy of such a 
chain is given by U e /k B T ~ q 2 (fL/A) 2 l B /y/(r 2 }- Elec- 
trostatic interactions become relevant for Ue/ksT > 1 or 
chain lengths L exceeding 



_ ;l/3 



^(7 



A 2 



fWlB 



2/3 



and whose spatial extension is given by 

A 2 ^'3 



i = i 



2/3 



f 2 q 2 h 



(3) 



(4) 



However, this derivation breaks down for strongly inter- 
acting systems where the contour length per blob becomes 
smaller than the intrinsic persistence length. In this case, 
a similar argument can be made for a WLC with L < Iq 
and (r 2 ) = L? yielding 



A 



q 2 Ph 



(5) 



Both definitions match for l g = lo, hiding a subtle crossover | 
behind a crudely renormalized system of units. Through- 
out the paper all quantities will be expressed using these 
natural units of contour length and spatial distance. On a 
scaling level, they become a function of only two dimen- 
sionless parameters: the reduced chain or segment length 



X = L/l g 

and the reduced screening length 



(6) 



(7) 

In the model under consideration the electrostatic in- 
teractions are purely repulsive. Therefore the chains are 



R. Everaers, A. Milchev and V. Yamakov: The Electrostatic Persistence Length of Polymers beyond the OSF Limit 3 



10 J 



10' 



GPVB 



BJ 



10 2 
10° 
10" z 







KMBJ 






^^^^^^.^ r a m 






QSFKK 













KK L 



10 



-2 



10° 10 : 



10* 



Fig. 1. Conformations of DHWLC beyond the OSF limit as a 
function of reduced chain length and reduced screening length. 
Electrostatic interactions are relevant in the upper right part 
of the map which is limited by solid black lines. The dashed 
and dotted lines correspond to crossover lines in the OSFKK 
and KMBJ scaling theories respectively. We have also included 
the scaling predictions Eqs. (^) to (j^j for the chain radii. For 
details see the main text. 



always extended relative to the neutral case. For kIb -C 
1 the details of the process are quite involved and dif- 
ficult to treat from first principles. The approximation 
schemes used are often based on mechanical analogies such 
as stretching, bending and swelling due to short-range 
excluded volume interactions and have been reviewed in 
Ref. |l|,§| • We have summarized the various predictions 
in terms of a schematic map of the XY parameter space 
(Figure [j]). Before we discuss the controversial parts, we 
first present those aspects which seem well understood: 

— The Debye-Huckel interaction is irrelevant inside the 
electrostatic blob, i.e. for weakly charged chains 



For sufficiently long chains, the Debye-Huckel interac- 
tion becomes effectively short-ranged, leading to self- 
avoiding walk (SAW) behavior 



(r 2 )/e 



X 1 



for < Y < oo and X -> oo (10) 



where v ss 3/5 is the usual Flory exponent. 
For strong screening with q 2 ls < k^ 1 < A the Debye- 
Huckel interaction reduces to an ordinary excluded 
volume potential with a second virial coefficient v ~ 
q 2 lsK~ 2 between charges. Using a conventional Flory 
argument to balance the two-body repulsion v(fL/A) 2 /R F 
with the entropic elasticity of a Gaussian chain R f /(LIq), 
one obtains: 



(r 2 )/e 



\ y 4 / 5 x 6 / 5 



for X < Y- 4 
for X > Y- 4 



for Y < 1 



(11) 



In Figure [l] the corresponding line, beyond which the 
short range excluded volume interaction becomes rel- 
evant, is marked as "F". 

The controversial parts of the phase diagram concern 
the crossover from the blob pole to the self-avoiding walk 
regime. The problem is often treated in analogy to a simple 
WLC. With an Onsager virial coefficient v ~ l 2 d between 
rigid segments of length l p and diameter d, excluded vol- 
ume effects become relevant beyond a "Flory length" lp = 
ll/d 2 leading to a Flory radius of R% ~ d 2 ' 5 f p /5 {L/l p f' 5 . 
In the case of the blob chain, the diameter of the elec- 
trostatically excluded volume is given by d ~ k" 1 0, 
p0{ . However, there is disagreement with respect to the 
K-dependence of the electrostatic persistence length l e . 

— Variational approaches such as the theory of Barrat 

d ~ As a 
_1 so that there is a direct 



and Joanny (BJ) often predict l e 



consequence, lp = l p - 
crossover from the stiff blob chain to a SAW regime 
when the contour length £ y of the blob chain reaches 
the screening length Using dimensionless units 

this corresponds to X = N/g = Zf/C = W£ = 1/( K £) = 
Y (the dotted line in Figure [I] marked "BJ"). The re- 
sult 



X for X < 1. 



(8) 



For strongly charged chains, this regime does not exist. 
In the absence of screening, when the monomers inter- 
act via an infinite range Coulomb potential, the chains 
are stretched into a "blob pole": 



(r 2 )/e 



X" 



for I < X < oo and Y — > oo (9) 



In Figure [l] we have marked the line dividing these two 
regimes as "GPVB" after the initials of the authors of 
Ref. H where the notion of the electrostatic blob was 
introduced. 



(r 2 ) /£ 2 ~ Y 4/5 X 6/5 for X > Y and Y > 1 (12) 



mar [H, 
Most t 



is identical to Eq. (|1 1|) . On a scaling level, the predic- 
tions of the BJ theory coincide with those of the ex- 
cluded volume theories of Katchalsky fl2] and Muthuku- 

nj[n). 

reories favour the relation l e ~ £ _1 k -2 first 
obtained by Khokhlov and Khachaturian (KK). KK 
argued that the OSF result Eq. (||) should also apply 
to a stretched chain of blobs with line charge density 
fql g /{M) so that Z e /£ = I/«) 2 = V 2 . Since l e » d, 
the resulting phase diagram is considerably more com- 
plicated. The blob chain starts to bend for reduced 
segment lengths X exceeding Z p /£ = Z e /£ = V 2 (the 
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"KKi" line in Figure |1|), while excluded volume ef- 
fects become relevant beyond If/£, = V( K £) 4 = Y A 
("KK 2 "). 



Y 2 X 



for Y 2 < X < Y 4 , ^ , 
4 for Y > 1 



y6/5 X 6/5 f orI> y 



(13) 

Finally, and in contrast to the KMBJ theory, the OS- 
FKK approach implies another crossover ( "KK3" ) within 
the SAW regime at Y = 1 from Eq. (£|) to Eq. @). 



3 Implications for Data Production and 
Analysis 

In general, scaling theories make two kinds of predictions: 
(i) about the existence of characteristic length scales or 
crossover lines in conformation space and (ii) about the 
asymptotic behavior of observables in the areas between 
these crossover lines. In principle, attempts at refutation 
can aim at either type of prediction. However, in the present 
case the identification of asymptotic exponents turns out 
to be particularly difficult. Apart from numerical prefac- 
tors and logarithmic corrections [] one is faced with four 
problems: 

— Although the various regimes predicted by the KMBJ 
and the OSFKK theory are characterized by different 
combinations of powers of X and Y , the exponents are 
often similar and the absolute differences between the 
predicted chain extensions relatively small. 

— At least on a scaling level all crossover lines meet at 
X = Y = 1 for chain and screening length of the order 
of the diameter of the polyelectrolyte blob. In the case 
of the OSFKK theory some of the predicted regimes 
are extremely narrow in the sense that chain lengths 
of X = 10 4 blobs are required for a width of one order 
of magnitude in Indirection. This validity range would 
seem to be the absolute minimum for identifying power 
law behavior. 



The crossover lines are neither parallel to each other 
nor to the "natural" X and Y directions of variations 
of chain and screening length respectively. 
In particular, results will be influenced by the finite 
total length L to t of polyelectrolyte chains studied in 
experiments or simulations. The importance of these 
effects varies with the ratio of the screening length 
Y and the contour length X of the blob chain. As 
a consequence, they risk to mask the asymptotic In- 
dependence of observables such as the electrostatic per- 
sistence length, if they are evaluated for chains with 
fixed Ltot- 

Quite obviously, the discrimination between the two scal- 
ing pictures requires the investigation of chains whose ef- 
fective length X is as large as possible. In addition one 
should rely on those observables and data representations 
which are most sensitive to the differences between the 
theories and least sensitive to the omitted constants, cor- 
rections and crossovers. In the following we discuss the 
analysis of data for internal distances and for the tangent 
correlation function. In particular, we will argue that it is 
relatively easy to discard the KMBJ picture using simple 
scaling plots, while the verification of some of the predic- 
tions by the OSFKK theory requires astronomical chain 
lengths. 

Consider first the scaling predictions Eqs. (g) to ( p^ ) 
for the mean square internal distances at reduced screen- 
ing lengths Y > 1. Describing the GPVB crossover to 
the blob pole within the chain- under-tension model [g,[32]], 
Eqs. (g) and (g) can be combined as (r 2 ) /£ 2 = X + X 2 . 
Taking this into account, the KMBJ theory suggests that 
all data points should collapse when plotted in the follow- 
ing manner as a function of the KMBJ persistence length: 



For example, in the absence of screening the long-range 
Coulomb interactions along the blob pole create a tension 
which grows logarithmically with the chain length TV. As a 
consequence, the blob diameter is reduced and the chains grow 
with R g ~ TV log 1//3 (TV). In particular, a TV monomer seg- 
ment of a longer chain will always be more extended than a 
TV monomer chain whereby the deformation is strongest for 
segments located near the center of the longer polymer [g|3l]] . 
In the presence of screening, this effect leads on the one hand 
to an increase of the contour length of the blob chain. On 
the other hand, the correspondingly reduced line charge den- 
sity results in a reduction of the OSFKK persistence length. 
Similarly, but neglecting the stretching, there are logarithmic 
corrections to the electrostatically excluded volume around the 
blob chain [Q. Moreover, although fairly robust, the Flory ar- 
gument used to estimate the excluded volume effects is far from 
exact. In addition to the aforementioned crossovers, a complete 
theory will have to account for all of these effects. 



(r 2 (Z = X/Y))/e-X 
X 2 



1 forZ<l 

Z- 4 / 5 for Z > 1 ti4) 



In contrast, the OSFKK theory predicts data collapse, if 
the segment lengths are rescaled with the OSFKK persis- 
tence length Y 2 , and a breakdown of scaling for segment 
lengths approaching the Flory length Y 4 : 



(r 2 (Z = X/Y 2 ))/e - X 
X 2 



1 for Z < 1 

Z- 1 for 1 < Z < Y 2 

y-2/5 Z -4/5 for Y 2 <Z 

(15) 

The predictions of the two scaling theories differ most 
strongly for chain radii along the Y — X 1 / 2 KKi line for 
segment lengths equal to the OSFKK persistence length. 
In the KMBJ theory, this line is already deep in the SAW 
regime. Eqs. (|l2|) and dill) imply 



(r 2 (X,Y ^X 1 ' 2 )) 

xe 



X 3 / 5 (KMBJ ) 
X (OSFKK 



(16) 



Thus the ratio (r 2 ) OSFKK / (r 2 ) KMBJ = X 2 / 5 is fairly 
small and X — IO 5 / 2 blobs are required for this ratio to 
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GPVB 




X=N/g 



Fig. 2. Areas of the conformation diagram Fig. [i] for which 
there are experimental and simulation data available. In (a) 
the green lines indicate the parameter ranges investigated in 

previous numerical studies: Barrat and Boyer |j2| ( ), Sei- 

del J33[ ( ), Jonsson et al. [ p4| , |28[ ( ), Micka and Kre- 

mer |p7|] ( ). Experiments (shown in red) have access to 

longer chains, but the reduced screening length Y are typically 

smaller than ten: Reed et al. ^] ( ), Beer et al. E| ( ). 

The colored grids in (b) denote the parameter ranges covered 
by different sets of our MC simulations. Note that the predic- 
tions of the KMB J and the OSFKK theory differ most strongly 
along the KKi line and are identical outside of the the gray 
shaded area. 



become of order 10. For comparision, both theories predict 
full extension along the KMBJ line 



{r 2 (X,Y = X)) 



X (KMBJ ) 
X (OSFKK ) 



(17) 



and differ only by a factor of (r 2 ) OSFKK / (r 2 ) 
X 1 / 10 along the KK 2 line 



KMBJ 



(r 2 {X,Y = X 1 ' 4 )) 

x~e 



X 2 / 5 (KMBJ ) 
X 1 / 2 (OSFKK ) 



(18) 



where Z is the variable measuring the effective distance 
from the crossover line at Z = 1. In order to identify the 
asymptotic behavior one needs at least data covering the 
interval Z e [0.1, 10]. Since the validity range of Eq. (|i~fj| ) 
is limited by the KKi and KK3 lines (so that 1 < Y < 
X 1 ! 2 or A" 1 / 5 > Z > AT -1 / 5 ), this implies a minimum 
segment length of X = 10 5 blobs for establishing the KK2 
crossover. Similarly the KK3 crossover between Eqs. (O) 
and (|l^) at Y = 1 becomes relevant for chains of X = lCr° 
blobs! 

The mean-square internal distances and the tangent 
correlation function (TCF) obey a Green-Kubo like rela- 
tion: 



(bjy • b ) = ~ 



2 -, f 2 j2 



(r 2 (X)) 



(20) 



2dN 2 ' y " 2g 2 dX 2 £ 2 
For a WLC the TCF is simply given by 

(b(s)-b(O)) =& 2 exp(-s/y (21) 

so that the persistence length can be read off directly 
from a semi -log arithmic plot. Numerical studies of poly- 
electrolytes |27|j28|,p9[ have therefore often focused on this 
quantity in spite of two intrinsic problems: (i) the TCF is 
considerably more difficult to measure with the same rel- 
ative precision than internal distances and (ii) the TCF 
is particularly sensitive to finite chain length effects (a 
characteristic sign is the faster than exponential decay of 
the TCF on length scale approaching the chain length). 
In contrast to the case of ordinary SAWs j55|, nothing is 
known about the functional form of the corrections. In the 
following discussion we will focus on a third aspect: the 
sensitivity of the TCF to the neighborhood of crossover 
lines. 

On a scaling level, the behavior of the TCF can be 
obtained by applying Eq. © to Eqs. (|) to (]l|). For 
Y > 1 the KMBJ theory predicts 



Ub(z = x/Y).b(o)) 



e 2 



for Z < 1 
4 / 5 for Z > 1 



(22) 



Within the OSFKK theory, simple predictions can only 
be made for segment lengths below the persistence length 
and beyond the Flory length: 



-(h(Z = X/Y 2 )-b(0)) 



While segment lengths of the order of X = 10 3 blobs 
are thus sufficient to discriminate between the KMBJ and 

the OSFKK proposals for the electrostatic persistence length, H °wever, since both theories are based on the analogy 
the requirements for resolving the additional crossovers 
predicted by the OSFKK theory are much higher. Con- 
sider again the KK2 line where excluded volume effects 
are expected to become relevant for the undulating blob 

chain. Eq. ( |l3| ) can be rewritten in the form j^(b(X) ■ b(0)) ~ (24) 



1 for Z < 1 

y-2/5 Z -4/5 for Y 2 <Z 

(23) 
to a 

mechanical WLC, they are often associated with the much 
more detailed prediction 



(r 2 (Z = y-VSAT 1 / 5 )) 

Y 2 xe 



1 for Z < 1 
Z for Z > 1 



(19) 



exp(- 
exp(- 



b(0)> 

X/Y) 
X/Y 2 ) 



for A < Y 
for A < Y 2 



(KMBJ ) 
(OSFKK ) 
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for the functional form of the decay of the tangent cor- 
relations. Measuring this quantity for DHWLC therefore 
seems to be the most direct way of justifying or refuting 
this analogy and its exploitation. In particular, numerical 
work |^,^8|,^9| has concentrated on (i) establishing the 
existence of a single exponential decay of the TCF over 
a certain range of length scales and (ii) extracting the 
K-dependence of the measured decay length. In the fol- 
lowing we will reexamine this approach by taking a closer 
look at Eqs. ( p2| ) and (|2^), since they contain additional 
crossovers neglected in Eq. (p4]). 

The situation should be uncritical for the GPVB crossover 
where the chain- under- tension models [g, 32 suggests that 
Eq. (|24|) remains valid for X < 1. In contrast, nothing is 
known in detail about the way the TCF crosses over to 
the slow power law decay characteristic for the SAW be- 
havior on large length scales. However, matching Eq. ( pi| ) 
(which only accounts for the local bendin g ri gidity) with 
the asymptotic behavior in Eqs. ( p2|) and ( |23| ) shows that 
the tangent-correlation function is much more sensitive 
to excluded volume effects than the chain radii. This is 
most obvious for the OSFKK theory where the two limits 
match close to the OSFKK persistence length X = Y 2 
instead of the Flory length X = Y 4 . While the scaling 
of the TCF with the OSFKK persistence length should 
start to break down around X/Y 2 w 1, one can neverthe- 
less expect Eq. ( pi| ) to hold up to this point. In the case 
of the KMBJ theory the situation is quite different, since 
the persistence length and the Flory length coincide. As 
a consequence, Eq. (^J) effectively breaks down as soon 
as the tangent-correlation function starts to deviate from 
one. On the other hand, in the absence of other relevant 
length scales the TCF should scale with the KMBJ per- 
sistence length for arbitrary segment length! 

In our opinion, these arguments shed some doubts on 
attempts to identify the electrostatic presistence length 
which are based too closely on Eq. (pi|). Scaling plots 
testing Eqs. ( |22| ) and ( p3|) may offer a simpler and safer 
alternative. 



case of a stretched chain with large excess charge. One 
MC step consists of N attempted rotations at ran- 
dom positions along the chain. Chain conformations 
are stored at intervals of 8-32 MC steps. For each pa- 
rameter set we simulate 8 independent Markov chains 
in parallel. We typically store 8 x 60 conformations 
representing a total of 1.5 x 10 7 attempted rotations 
for our longest chains. 

Instead of the slower procedure of Stellman and Gansj38| 
which corrects the accumulating numeric error in off- 
lattice implementations of the pivot algorithm, we reg- 
ularly reconstruct the chains with the correct bond 
length. 

— For calculating the long-range electrostatic interactions 
we use a direct summation whereby the energy of the 
system is obtained by direct counting of all the pair 
energies of the beads. This method is still efficient for 
macromolecules of up to few thousand monomers. The 
DH potential is tabulated in two arrays for short and 
long distances respectively. 

— For better efficiency starting configurations of the chains 
are created by means of the configurational biased [|39| 
MC method although one should keep in mind that 
due to the long range interactions the first part of the 
newly grown chain does not experience the cumulative 
field of the rest of the chain and a number of rotational 
moves are still needed before the chain is well equili- 
brated. Measurements are performed and conforma- 
tions stored only after the chain end-to-end distances 
are well equilibrated. 

— Since statistics is gathered both with respect to chain 
conformations as well as to different Bjerrum length 
Ib and screening length k , we use a simple paral- 
lelization where different processors of a CrayT3E su- 
percomputer perform independent simulation of single 
chains. The total CPU time used for this project is of 
the order of 1.510 5 single processor hours. 

The simulation parameters translate into our blob units 



4 Simulation Model, Method and Parameters 

As already mentioned in section 0, we model the poly- 
mers as freely jointed chains (FJC) with unit charges q = 
1 at each joint. Lengths were measured in units of the 
bond length b. We varied the the Bjerrum length Ib = 
4 2 , 1, 1/4 2 , 1/16 2 , 1/100 2 6 and the screening length K" 1 = 
1, 2, 4, 8, 16, 32, 64, 128 b. As in our previous study on poly- 
ampholytes j$6|, the chains have a length of up to W = 
4096 monomers. 

— Since we study the conformations of isolated chains, we 
employ the efficient technique of pivot rotations due to 
Sokal et al.|37,34|. We use two types of pivot moves: 
Either we rotate the part between the free end of the 
chain and a randomly selected monomer around an 
axis, defined by the bond between this monomer and 
its nearest neighbor; or we rotate a segment between 
two randomly selected monomers around an axis join- 
ing them. The latter provides better efficiency in the 



2/3 



(s > 1) 



1/3 



(25) 
(26) 



and 



T- (-9 < 1) 



(27) 
(28) 



where we now use the number g of monomers per blob 
instead of the corresponding contour length l g = bg. 

Fig. (|2|) shows where our own data are located within 
the AF-conformation space. The effective chain and screen- 
ing lengths studied cover a range of seven and five orders 
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Fig. 3. Comparison of measured internal distances {r 2 (X, Y))/^ 2 to the predictions of the KMBJ ((a) and (c)) and OSFKK 
((b) and (d)) scaling theories. In the top row we show log-log-log representations where all distances are normalized to the 
undisturbed random walk. The colored areas were generated by interpolation between the results of all simulations for a given 
coupling strength. The supporting grids and the crossover lines show the two sets of scaling predictions Eqs. ([]) to (O) as an 
extensions of Fig. (Q) to three dimensions. In the plots of the bottom row the colors indicate the ratios (r 2 (X, Y)) / (r 2 ^~ 

and (r 2 (X,Y))/ (r 2 ) KMB J respectively. 



' OSFKK 



of magnitude respectively. The reduced mean square inter- 
nal distances vary over ten orders of magnitude. Along the 
KKi line our data extend on a logarithmic scale about a 
factor of two further into the asymptotic regime than pre- 
vious studies. While this allows us to discriminate between 
the KMBJ and the OSFKK predictions for the electro- 
static presistence length, our chains are still too short to 
resolve the different RW and SAW regimes predicted by 
the OSFKK theory. 

Note, that only by studying strongly stretched chains 
we are able to push the effective chain length X close to 
10 5 and that our unified description of strongly and weakly 
charged flexible polyelectrolytes needs to be confirmed by 
the data analysis. To facilitate the comparison, all figures 
make use of the same color code to indicate data obtained 
for a particular coupling strength Is/b ranging from blue 
for g = 10000 2 / 3 « 470 over different shades of violet for 
g = 256 2/3 w 40 and g = 16 2/3 w 6.4 to red for g = 1 and 
orange for g = 1/16. The first three systems can safely be 
regarded as Gaussian chains, while the last two are at and 
beyond the crossover to the strong stretching regime. 



5 Results 

In the data analysis we mainly concentrate on identify- 
ing the scaling behavior: (i) Do different data sets overlap 
when rescaled according to our extension Eq. (g) of the 
definition of the electrostatic blob? (ii) How do the results 
of our simulations compare to the predictions of the KMBJ 
and OSFKK scaling theories? In terms of observables we 
start by presenting data for internal distances averaged 
along our chains of total length N — 4096. In the second 
part, we discuss results for the tangent correlation func- 
tion. While the TCF was also averaged along the chain, 
we only take into account distances up to half the chain 
length in order to reduce finite-chain length effects |p5[ . 
Except for the most weakly charged system, the chain ex- 
tensions are much larger than the screening length, so that 
we do not expect finite chain length effects to be very im- 
portant. At the end, we briefly present results for shorter 
chain lengths. 

Figs. ||(a) and (b) are three-dimensional log-log-log 
plots giving an overview of all data. We show reduced in- 
ternal distances (r 2 (X,Y)) / ((N/g)^ 2 ) normalized to the 



8 R. Everaers, A. Milchev and V. Yamakov: The Electrostatic Persistence Length of Polymers beyond the OSF Limit 



10* 
ID* 
10 2 

1 











BJ: Y=X . 




KK ± : Y=X 1 ^_ 








KK 2 : Y=X 1 ^ 




10" 2 3 


x 10 2 10 J 



Fig. 4. Extension of chain segments along the BJ, KKi, and 
KK2 crossover lines in comparison to the predictions of the 
OSFKK (solid line) and the KMBJ (dashed line) scaling theo- 
ries (see Eqs. (0), (^6|), and (|l8|); note that we have not used 
additional prefactors for this comparison) . Results for different 
coupling constants are shifted by factors of vTOOO- 
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size of the undisturbed random walk as a function of the 
reduced chain and screening lengths Eqs. ([|) and (0) re- 
spectively. Results for different coupling constants are com- 
bined into colored surfaces, while the supporting grids 
show the two sets of scaling predictions Eqs. to ( [l3| ) 
as an extensions of Fig. |l| to three dimensions. 

The complementary Figs. ||(c) and (d) show the ratios 
(r 2 (X 7 Y))/(r 2 (X,Y))KMBj and (r 2 (X, Y))/(r 2 (X, Y)) OSFKI 
of the interpolated simulation results to the scaling predic- 
tions in a color coding where green, red and blue indicate 
agreement, under- and overestimation by a factor of three 
or more respectively. The advantage of this representation 
is the localization of deviations in our schematic map of 
the parameter space. 

Qualitatively, the interpretation of Fig. |j| seems clear. 
There are neither indications for a failure of the blob scal- 
ing nor for significant deviations from the predictions of 
the OSFKK theory. (We emphasize again that we have 
neglected all numerical prefactors and that Eqs. (|J) to 
(O) treat crossovers in the crudest manner). In particular, 
there is no evidence that the chains start to bend on length 
scales comparable to the screening length as predicted by 
KMBJ . In the relevant part of conformation space, the 
KMBJ theory systematically underestimates the observed 
chain extensions. 

Nevertheless, Fig. || could be misleading, since the re- 
jection of the KMBJ theory is mainly based on data falling 
into the strong-stretching regime, while the theory is meant 
to apply to weakly stretched Gaussian chains. Thus so 
far our conclusions rest on the assumption that the ex- 
tension Eq. (||) of the blob scaling to strongly charged 
chains can be used to extrapolate the behavior of weakly 
charged systems to segment lengths inaccessible by sim- 
ulation. How well this assumption is fulfilled is hard to 
judge from Fig. [|. Definite conclusions require a more de- 
tailed analysis. 

Fig. ^ presents chain radii measured along the BJ, 
KKi and KK2 crossover lines. The first point to note is 
that in all three cases we observe almost perfect scaling of 




X/Y 1 X/Y 1 



Fig. 5. Crossover scaling for internal distances versus segment 
length. In the top (bottom) row segment lengths X are nor- 
malized to the KMBJ (OSFKK ) persistence length Y (Y 2 ) 
respectively. Figures (a) and (c) on the left-hand side show 
internal distances {r 2 (X, Y)) normalized to the mean-square 
extension £ 2 X 2 of the blob pole. The grid is the same as in 
Fig. |^(b) and shows const — X and const — Y lines. Figures 
(b) and (d) on the right-hand side are inspired by the chain- 
under-tension model for Gaussian chains Eqs. ( p"i| ) and (|l5|). 
Only const — X lines are shown. Data points falling into the 
range 1 < Y < X < Y 2 are marked using the color code in- 
dicating the coupling strength. Results for different coupling 
constants are shifted by factors of ten. 



data obtained for different coupling constants. Eqs. ( |16| ) 
to (|l^) predict simple crossovers at the blob size around 
X = 1. In agreement with both scaling pictures, we ob- 
serve stretched blob chains along the BJ line. The most 
important set of data are the radii measured along the 
KKi line for chains with a contour length X — Y 2 equal 
to the OSFKK persistence length. In agreement with the 
OSFKK theory we find a simple crossover around X = 1 
to stretched blob chains. Contrary to the predictions of the 
KMBJ theory the radii are essentially identical to those 
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Fig. 6. Scaled tangent correlation functions. Results for differ- 
ent coupling constants are shifted by factors of ten. The colored 
lines mark the results of our fits to a simple exponential decay 
in the range Y < X < Y 2 . (a) log- log plot using KMBJ scaling, 
(b) log-log plot using OSFKK scaling, (c) semi-log plot using 
OSFKK scaling. 



observed along the BJ line and do not show SAW behav- 
ior. In particular, the asymptotic slope predicted by the 
OSFKK theory is already observable for weakly charged 
chains to which the KMBJ theory can be applied directly. 
The last set of data is taken along the KK2 line which 
marks the onset of excluded volume effects in the OSFKK 
theory. Here our results are consistent with the predictions 
of both theories. This observation is in agreement with the 
estimate of a minimum segment length of X = 10 10 blobs 
for the difference to become relevant (see Eq. (p!s|)). 

The screening length dependence of the effective bend- 
ing rigidity of the blob chain can also be determined from 
crossover scaling of internal distances normalized to the 
size of the stretched blob pole (Fig. |^). The disadvantage 
of plotting (r 2 {X,Y))/(£_ 2 X 2 ) directly (Figs. | (a) and (c)) 
is the occurence of a 1/X divergence of results for segment 
lengths smaller than the blob size. Correcting for this in 
the manner su gges ted by the chain-under-tension model 
(Eqs. ( [Til ) and (|15|)) as in Figs. || (b) and (d) largely elim- 
inates effects due to the GPVB crossover, but introduces 
some artifacts for small N where (r 2 (N = 1)) — b 2 e for a 
FJC. In agreement with our previous results and indepen- 
dently of the coupling constant we observe extremely poor 
scaling when the data are plotted as a function of the ra- 
tio X/Y of chain length over KMBJ persistence length. In 
constrast, the data superimpose considerably better, if the 
OSFKK scaling is used as in the corresponding Figs. |^(c) 
and (d). In particular, Fig. |^ eliminates the possibility 
of an electrostatic persistence length scaling like n" 1 but 
with an unusually large prefactor. 

Similar conclusions can be drawn from an analysis of 
the tangent-correlation function (Fig. [3]). In fact, our dis- 
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Fig. 7. Location of apparent electrostatic persistence lengths 



le 



in our schematic map of the XY parameter space. We 



show results for the reduced crossover distance X c d = l c d/£, (a) 
and the reduced orientational correlation length X oc — Z oc £ (b) 
for chains of length Ntot = 4096 (solid lines) and N to t = 256 
(dashed lines) . In the first case, we compare (r 2 ) {X,Y) to the 
size of the stretched blob pole Eq. fil and define X c d = l c d/£ 
implicitely via (r 2 ^ (X c d,Y) = ^ 2 X cd /3. The results for X oc 
presented in (b) are decay lengths extracted from fits of TCFs 
to simple exponentials. For the fits we used data from the seg- 
ment length interval X < Y, \Ntotjg- However, for Y < 10 the 
decay of the TCFs ceases to be well described by a simple ex- 
ponential (see also Fig. ^). Being strongly depend on the data 
range selected for the fits (data not shown, the values presented 
in (b) for Y < 10 thus have to be taken with a grain of salt. We 
note that the results for datasets with different coupling con- 
stant scale quite well and that there is good qualitative agree- 
ment between the two methods. The results nicely follow the 
OSFKK prediction for reduced screening lengths Y > 10, but 
are strongly influenced by excluded volume effects for smaller 
values of Y. In particular, the extracted persistence lengths 
systematically exceed the OSFKK estimate. With respect to fi- 
nite chain length effects the first method turns out to be more 
robust than the second. 



cussion in section || shows that Figs. || (b,d) and the log- log 
plots in Figs. || (a,b) are directly comparable. Fig. ^ (c) 
shows the same data in the semi-logarithmic representa- 
tion commonly used to identify a simple exponential de- 
cay of the correlation function. Clearly, the OSFKK scal- 
ing does not work perfectly up to the OSFKK persistence 
length, but, at least qualitatively, we observe the expected 
slow-down of the decay of the correlations. 

The discrepancy between our conclusion, l e ~ k~ v 
with y = 2, and the results of previous numerical and 
experimental investigations, y <C 2, can be traced back to 
the definition of the electrostatic persistence length. So far 
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we have used an indirect method (scaling plots of segment 
radii and TCFs). More direct methods usually proceed by 
(i) defining an apparant electrostatic persistence length 
le,app, (h) calculating l e , a pp for numerical or experimen- 
tal data, (hi) plotting l e> wm>/£, as a function of the reduced 
screening length Y (Fig. Q), and (tentatively) (iv) extract- 
ing effective values for y. Suitable definitions for l e . a pp were 
recently reviewed by Ullner and Woodward . Fig. [[] lo- 
cates our results for the "orientational correlation length" 
l oc and the "crossover distance" l cc i in our map of con- 
formation space. The first length is defined as the decay 
length of a simple exponential fitted to the TCF while the 
second tries to identify the crossover from the blob pole 
Eq. (H) to an undulating blob chain Eqs. ( |l^ ) or (|l3|). 

Clearly, the extracted length scales can only be identi- 
fied with the electrostatic persistence length l e as long as l e 
is well separated from other relevant length scales, i.e. for 
sufficiently large chain and screening lengths. Fig. [7] shows 
that our results closely follow the KKi line for N = 4096 
and Y > 10. However, the violation of either of the two 
conditions leads to deviations. Relative to the OSFKK 
prediction the extracted persistence lengths 

increase for small reduced screening lengths due to ex- 
cluded volume effects and 

decrease in the opposite limit due to finite chain length 
effects which are particularly strong for the TCF and 
quantities related to it. 

Depending on the definition of l e ,app and the range of chain 
and screening lengths studied, the combination of these 
two effects can lead to the observation of effective expo- 



nent L 



v which are much smaller than 



y = 1. 

However, since the weak n dependence of l e ,app is an ar- 
tifact of the definition of the quantity, there seems to be 
no contradition to the OSFKK theory. Attempts along 
these lines |27|,|28|,^9| therefore risk to create more con- 
fusion than insight as long as l e , a pp is not defined within 
a theoretical framework which explicitly accounts for ex- 
cluded volume effects 



6 Discussion 



indeed possible to understand DHWLC by considering a 
hierarchy of effects due to interactions between different 
classes of monomer pairs: 

Stretching due to the (effectively unscreened) Coulomb 
repulsion between neighboring monomers into a chain 
of blobs which has a finite 

Bending rigidity due to the screening of interactions be- 
tween monomers with a distance larger than g/(n£,) 
along the chain. As a consequence, the blob chain re- 
mains straight up the electrostatic persistence length 
l e = n~ 2 /t;. Beyond l e the chain behaves like a random 
walk, before 

Swelling due the electrostatically excluded volume between 
chain segments with a distance larger than l e becomes 
relevant beyond the Flory length lp = k~ 4 /£ 3 

An interesting side result of our work is the extension 
of the polyelectrolyte blob scaling to the case of strongly 
interacting, almost fully stretched chains for which the 
OSF theory is known to work This extension was done 
in the logic of the OSFKK theory but is incompatible with 
the ansatz of BJ. Its success provides strong evidence for 
the irrelevance of longitudinal and transverse fluctuations 
within the blob chain jl!| and proves the KK idea almost 
by itself. 

Clearly, scaling arguments cannot do justice to the full 
complexity of the problem. Omitting all numerical prefac- 
tors, the ubiquituous logarithmic corrections, finite chain 
length effects and, in our opinion most importantly, a re- 
fined description of the crossovers between narrow neigh- 
boring regimes, they cannot hope (and should not be ex- 
pected) to describe numerical or experimental data in de- 
tail. Quite obviously, these features call for a quantita- 
tive explanation. While our numerical results can serve as 
benchmarks for the development of theories, it is a sober- 
ing thought that the simplest model of a single, isolated 
polyelectrolyte chain is still unsolved. Compared to the 
much better understood neutral polymers, the OSFKK 
theory represents the equivalent of the standard Flory ar- 
gument for the excluded volume effect. Nevertheless, we 
believe to have shown that the OSFKK theory provides 
the indispensible "big picture" needed for the design and 
analysis of experiments and computer simulations. 



In this paper we have combined a scaling analysis of the 
conformational properties of intrinsically flexible polyelec- 

trolytes with Debye-Hiickel interactions with extensive Monte7 Acknowledgements 

Carlo simulations of isolated chains. Our study was fo- 
cused on the controversial case of polyelectrolytes beyond 
the OSF limit, i.e. to the case where the electrostatic 
screening length k _1 exceeds the bare persistence length 
of the polymers in the absence of electrostatic interactions. 

Our main result is the refutation of theories ||[l3| pre- 
dicting an electrostatic persistence length scaling as 
In contrast, we have observed no significant deviations 
from the scenario proposed by Khokhlov and Khacha- 
turian J^J who combined the idea by de Gennes et al. |^| of 
a stretched chain of polyelectrolyte blobs with the Odijk- 
Skolnick-Fixman theory of the electrostatic persistence 
length HH and the electrostatically excluded volume || 
|30|] between chain segments. Our results suggest that it is 
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